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Abstract 
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, Probability for a primordial black hole to invade the Kuiper belt was calculated. We 

pP showed that primordial black holes of certain masses can significantly change asteroids' 

orbits. These events may result in disasters, local for our solar system and global for the 



Earth (like the Tunguska meteorite). We also estimated how often such events occur. 



1 Introduction 



Nowadays the asteroid danger is one of the most important scientific issues. Geophysical 
data indicate that in the past the Earth was dangerously attacked by asteroids about 
every 190 million years rather than continuously. These repeated global events suggest 
that the probability for asteroids, or meteorites, to shoot is not constant, but depends on 
external (with respect to the Solar system) factors (see pQ). 

At the moment the nature of dark matter is not clear. One of the candidates for dark 
matter particle is a primordial black hole (PBH). In case the entire dark matter consists of 
PBHs, it is possible to set not only the lower (Hawking) bound for their masses, but also 
the upper one. The latter comes from the data on gravitational microlensing (for example, 
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■ see [2]). However, the upper bound is quite uncertain and in different estimations lies in 

the range 10 26 -r 10 28 g, i.e. close to the Earth and Moon masses. 

In this paper we obtained the probability for a PBH travelling through the Solar 
system to collide (interact) with an asteroid. This process can make collision of the Earth 
^ . and some heavenly body more probable. 

To date, we are aware of several asteroids whose orbits are threateningly close to the 
Earth. The asteroids that are sufficiently big and massive are most dangerous to our 
planet. For example, the Tunguska meteorite must have been around 50 m in diameter. 
As is well-known, there are three prominent asteroid belts in the Solar system: 

1. Mass of the main asteroid belt (the Phaeton belt) is about 10 24 g ~ 10~ 3 M^. 
Asteroids populating the belt are 2.5 — 3 astronomical units (AU) away from the Sun. 

2. Mass of the Kuiper belt is around 10 27 g ~ M®. Its asteroids are at the distance of 
40—60 AU from the Sun. 

3. The Oort cloud is located on the Solar system border (more than 100 AU), and its 
mass estimations are rather vague. 

As is known, probability to collide with a heavenly body is in direct proportion with its 
surface area. Total mass of the asteroids is comparable with mass of the Earth, but since 
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there are so many of them, their overall surface area appears to be almost equal to the value 
for the Sun. Besides, there is one more factor which increases the collision probability. An 
asteroid mass is considerably smaller than that of the Sun. This implies that collision with 
a PBH should not be necessarily direct to change the asteroid's trajectory. In order to 
accomplish this, the PBH should only be gravitationally scattered on the asteroid, i.e. it 
should pass by the asteroid at some distance. This phenomenon increases the probability 
for a PBH and an asteroid to interact as much as for several orders of magnitude. 

Thus, PBHs interact much more probably with asteroids than with the planets or even 
the Sun. 



2 Probability to find a PBH in the Solar system 

According to our hypothesis we suppose that PBHs are responsible for almost entire 
mass of our Galaxy. Simply put, we set the same mass m p bh for every PBH (in a more 
complicated case one can just introduce a distribution function of the PBH mass m^). 

The PBH concentration in the Galaxy is n p bh = Pgai/rripbh, where p ga i ~ 10~ g/cm 3 
is the average mass density in the Galaxy. 

Therefore, the frequency of the Solar system - PBH collisions is 

v = ttRI v p ga i/m pbh , (1) 

where R a is the average radius of an asteroid orbit, and v is the average asteroid velocity 
in the Galaxy (at infinity with respect to the Solar system). 



3 Probability for a PBH to collide with an aster- 
oid 

Let the total mass of the asteroids be M as t. Density of an asteroid is usually in the 
following range: p as t = 2-^7 g/cm 3 . Asteroid sizes cover a large interval, with most 
of the rocks having less than 1 km in diameter and smaller size asteroids being more 
abundant. Since small asteroids do not result in significant impact, they are not interesting 
to deal with. So it makes sense to consider asteroids bigger than some critical size, e.g. 
D as t ~ 100m. One can estimate their number as 

N = last ■ Mast (2) 
Past • Dl st 1 ' 

where j a st is a mass fraction of the asteroids bigger than D as t in diameter. 

The frequency v of the PBH - asteroid interaction is obtained by multiplying the 
frequency vq by the probability v\ for the PBH to interact with a Solar system asteroid, 
which is 

v\ = Sast/(nRl) , Sast = N ast ■ irh 2 , (3) 

where h is the maximum impact parameter when the PBH changes significantly the as- 
teroid trajectory. 

Hence, using eq. ([1]) we obtain: 

T^last PgalM ast vh 2 

v = o . (4) 

m pbh p ast Dl st 

The quantity h is to be calculated in the next section. 
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4 Abrupt change of asteroid orbits due to scat- 
tering on PBHs 

With the PBH - asteroid interaction being the classical two-body problem, it is more 
convenient to solve it in the center of mass frame. 

In this reference frame, it is convenient to introduce the following notation^: 
m = m ast rripbh/ ( m ast + m phh ) is the reduced mass, 
P = m v is the reduced mass momentum in the center of mass frame, 
L = mvh is the angular momentum of the PBH - asteroid system, 
tq is the minimal distance between the PBH and the asteroid, 

Uq = Gm p bh m ast/ r o is an absolute value of the gravitational interaction energy in the 
orbit's perihelion, 

E = \mv 2 is the total energy of the system, 

A = Uq/ (2E) = G (rripbh + ~m ast ) / [tq v 2 ) is a dimensionless scattering parameter, 
B = ro/h is the ratio of the minimal distance to the impact parameter. 

Subtracting the zeroth order n and taking into account relations between the orbit 
parameters, we obtain the angle of the momentum deflection after the PBH has been 
scattered by the asteroid [3]: 



A(p = 2 arccos [1] — 2 arccos 
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As for B, the calculation yields: 
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In the center of mass frame, the change of the momentum is H: 

(AP) 2 = 2P 2 • [1 - cos(A^)] . 



Using the simple trigonometry and eqs. 

AP = P 2A 



5]) and ([7]), we obtain: 

2Gm(m pbh + m ast ) 



l + A 



2PA 



r v 



(5) 



(6) 



(7) 



(8) 



In the Galaxy the PBH average velocity is comparable to the star velocity disper- 
sion: v ~ 300 km/s. Therefore, the total energy E of the system is much more than the 
gravitational interaction energy: E » Uq, hence, A « 1, B ~ 1. 

The condition of the crucial alteration of the asteroid orbit is equality between the 
initial asteroid momentum, P as t = rn as tv as t, and the gained momentum AP within an 
order of magnitude: 

2PA 

aP ast = a m ast v ast = AP = - - — , (9) 



where, according to eq. ([7|), the parameter a = \/2(l — cos Aip) reflects how strongly the 
asteroid is affected. 

Taking into consideration that the asteroid velocity on a circular orbit of radius R is 
v as t = \J GMq I R and making use of eq. ([9]), we obtain the exact expression for h 2 , 
mass of asteroid influence having been taken care of: 



the 



2 = ^ Gm lbh R 
1 a 2 v 2 M 



a(m pbh + m ast )v ast 
2m pbh v 



(10) 



1 Here, G is the gravitation constant, M Q is mass of the Sun, M 9 is mass of the Earth. 

2 An absolute value of the momentum change is the same both for the asteroid and the PBH. 
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or 



''ast 



•800 km, (11) 



where R as t ~ 10 15 cm is the average distance to the Kuiper belt. 
Using this last but one expression, we obtain for A: 



a 



l GM Q I m pbh + m ast 
4v 2 R y m pbh 



« 1 (at m pbh > m ast ) . (12) 



This last expression validates the above-made assumption on A's being small for most 
asteroids of the Kuiper belt. 

Eqs. @ and (fTUj) yield the probability for an asteroid to be deflected by a PBH 
provided the latter has gotten into the Solar system: 

^ 4Gm 2 pbh j ast M ast ^ 2 ■ last ■ m 2 pbh ■ (100m) 3 

Substituting eq. (fTOj) to eq. (jH), we obtain the frequency of significant orbit deviations 
for big asteroids: 

_ 4 7T G J ast Pgal M ast R m phh 



a 2 Mq p as t Dl st V 



(14) 



5 Conclusion 

It is convenient to rewrite the final expression for the frequency of the events under 
consideration as follows: 

( lV- 7 ys- l \ (m pbh \ (M ast \ /t^IOOuA 3 / R \ 

v ~\^^)'Vm^)'Vm^) \— DaZr) \r^)' (15) 

The last three factors in formula (|15p are of the order of unity, so that the final expression 
only depends on the distribution and total mass of PBHs in the Galaxy. 

Estimating the fraction of asteroids hitting the Earth itself is quite difficult. This is 
due to the fact that, even if the asteroid does not hit some planet, very probably it will 
not leave the Solar system. Its orbit now crosses orbits of planets and other asteroids 
and becomes more dangerous. So sooner or later the asteroid itself or secondary asteroids 
(produced by the very first one) fall onto the planets. 

Thus, the probability for the scattered asteroid to hit the Earth can be roughly esti- 
mated as the ratio of the Earth surface area to the total surface area of the other planets. 
It yields approximately ~ 10~ 2 . Hence, it does not seem improbable that formula (|15p 
may give hit periods of about 190 million years, which is consistent with observational 
data. 

If the hypothesis analyzed in this paper is correct, modern methods aimed at averting 
the asteroid danger appear to be inefficient. This is related to the fact that their main 
idea is revealing big meteors and asteroids with dangerous orbits and, then, monitoring 
these bodies. However, if the main danger consists in abrupt changes of asteroidal or- 
bits (because of scattering on a PBH), revealing potentially dangerous bodies is hardly 
possible. 
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